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ABSTRACT 

A survey of CO(3-2) emission from optically visible oxygen-rich Mira variable stars within 500 pc 
of the sun was conducted. A molecular envelope was detected surrounding 36 of the 66 stars examined. 
Some of these stars have lower outflow velocities than any Miras previously detected in CO. The average 
terminal velocity of the ejected material was 7.0 kms~', about half the value found in Miras selected by 
infrared criteria. None of the stars with spectral types earlier than M 5.5 were detected. The terminal 
velocity increases as the temperature of the stellar photosphere decreases, as would be expected for a 
radiation driven wind. Mass loss rates for the detected objects were calculated, and it was found that 
there is no correlation between the infrared color of a Mira variable, and its mass loss rate. The mass 
loss rate is correlated with the far infrared luminosity, although a few stars appear to have extensive dust 
envelopes without any detectable molecular wind. A power-law relationship is found to hold between 
the mass loss rate and the terminal velocity of the ejected material. This relationship indicates that the 
dust envelope should be optically thick in the near infrared and visible regions of the spectrum when the 
outflow velocity is ^ 17 kms~'. At the low end of the range of outflow velocities seen, the dust drift 
velocity may be high enough to lead to the destruction of the grains via sputtering. Half of the stars 
which were detected were re-observed in the CO(4-3) transition. A comparison of the outflow velocities 
obtained from these observations with those obtained by other investigators at lower frequencies shows 
no evidence for gradual acceleration of the outer molecular envelope. 
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1. Introduction 



Over the last 25 years the new infrared surveys, particularly the Two-Micron Sky Survey (Neugebauer 
& Leighton 1969, hereafter TMSS), the Air Force Geophysical Laboratory survey (Price & Walker 1975, 
hereafter AFGL), and the IRAS Point Source Catalog (Neugebauer et al. 1984, hereafter PSC), have 
revealed the presence of tens of thousands of cool stars surrounded by dusty envelopes. After each of 
these surveys was completed, observers examined the most prominent of the newly detected objects in 
the radio regime, and were frequently able to detect molecular line emission from molecules such as OH, 
HCN and CO. From this confluence of infrared and radio data arose a compelling model for mass loss in 
highly evolved stars. According to this model, dust condenses in the outer envelope of a cool red giant 
star. This dust is then accelerated by radiation pressure to a velocity of ~ 15 kms~'. The dust in turn 
accelerates gaseous material by coUiding with it (see a review by Olofsson 1988). 

This paper presents the results of a survey of Mira variables with distances of 500 pc or less, which 
were selected without regard to their infrared characteristics. There are several reasons for carrying out 
such a survey. First, it is possible that the radiation pressure mechanism described above is not the only 
one driving mass loss from red giant stars. Mechanisms which do not require the formation of dust, such 
as radial pulsations (Wood 1979), momentum transfer from Alfven waves (Hartmann & MacGregor 1980) 
or from sound waves (Fusi-Pecci & Renzini 1975, Pijpers & Heam 1989) have been proposed for driving 
AGB winds. Finding a large, expanding molecular envelope surrounding a highly evolved star without 
a massive dust envelope would be strong evidence for an alternative mass loss mechanism. Secondly, 
because previous surveys tended to concentrate on the very strongest infrared sources, stars with small, 
but detectable, mass loss rates may be under-represented in their compilations of detections. The final 
reason for surveying nearby Miras is that they obey a period - luminosity relationship. This allows their 
distances to be estimated more accurately than many of the other types of stars which have been detected 
in the surveys of bright infrared objects. The derived values of many important physical parameters, 
including the absolute luminosity and the mass loss rate, are sensitive to errors in the estimated distance. 

Examining the CO(3-2) line, rather than a lower J transition, offers several advantages. The first 
is that higher J transitions have higher critical densities for colUsional excitation. For this reason the 
emission from molecular clouds is less extensive, and the spectra of Miras near the Galactic plane are less 
apt to be contaminated by galactic emission. Secondly, the antenna temperature of the 3-2 line will in 
general be higher than in the lower J lines (Morris 1980). While the rotational transitions above 3-2 may 
have even higher brightness temperatures, particularly for optically thin envelopes, the earth's atmosphere 
is much less transparent where these transitions are found, even on high mountain sites. The 3-2 line 
may therefore allow the most sensitive ground-based searches for cool circumstellar envelopes. Finally, 
for those objects which have akeady been detected in lower J surveys, the availability of spectra covering 
several CO transitions allows the models of molecular excitation to be better constrained, and variations 
of the CO envelope at different radii can be studied (van der Veen & Olofsson 1989). 

2. Observations 

The candidate objects for this survey where selected from the Revised Catalog of Spectra of Mira 
Variables of Types Me and Se (hereafter RCMV) (Keenan et al. 1974). The spectral types and periods 
of the stars were used to calculate an absolute visual magnitude, using the formulae derived by Cells 
(1980). This calculated absolute magnitude was used with the apparent magnitude, to calculate the star's 
distance. The star' s sky position and distance were then used to calculate the amount of absorption, using 
Parenago's formula with coefficients from Sharov (1964). Finally the amount of absorption was used to 
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derive a new distance estimate. This process of estimating the distance was iterated until it converged. 
Distances to a smaller set of stars were calculated in this way by Cells (1981), who estimated the mean 
error in derived distance values to be 19 %. However distances to many of these stars were later calculated 
by Jura and Kleinmann (1992) using infrared data, and their distance estimates differ from those of Celis 
by 50% on average. 79 of the stars in the RCMV were found to have a distance of 500 pc or less. 66 of 
these stars were observed. 

All observations presented here were made at the Caltech Submillimeter Observatory,* (CSO) on 
Mauna Kea, Hawaii. Data were taken during January, May and August of 1990, and during February 
through August of 1991. The i2ci6o(3-2) hne at 345.796 GHz was observed with an SIS receiver (ElUson 
& Miller 1987) having a double sideband temperature of about 200 K. A 1024 channel AOS spectrometer 
was used, which has a total bandwidth of 500 MHz and a velocity resolution of 0.9 kms~^. In most 
cases, scans were taken while position switching at 0.05 Hz between the source position and 2 off-source 
positions located 3 arc minutes away in azimuth. For weak objects at low galactic latitudes, the scans used 
4 off-source positions spaced symmetrically about the source, and 30 arc seconds away (1.5 beam widths), 
in order to remove any contamination from extended galactic emission. For those stars which appear in the 
SAO catalog, SAO coordinates and proper motions were used. For the remainder of the stars, coordinates 
were obtained from the IRAS PSC. Table 1 lists the names, coordinates and on-source integration times 
for the stars. For the stars which had not previously been detected in a rotational transition of CO, the 
radial velocity used for receiver tuning was obtained from the RCMV. If the star showed only weak CO 
emission (Tmb < 200 mK), it was re-observed on at least one additional night in order to confirm the 
detection. Two of the stars which were not detected, R Car and R Cen, never rise above 10 degrees 
from the horizon when observed at the CSO. Since telescope pointing can be problematic at very low 
elevations, the pointing was checked after these two stars were observed by observing the strong source 
IRC+10°216 as it passed through an elevation of 10 degrees. Both the CO(3-2) and CS(7-6) transitions 
were seen, indicating that errors in pointing were not responsible for these non-detections. 

In April - May of 1993, and January 1994, nineteen of the stars with strong CO(3-2) emission were 
examined again with a second SIS receiver (Walker et al. 1992), tuned to the CO(4-3) line at 461.041 
GHz, where it had a double sideband temperature of ~ 250 K. The system temperature was ~ 2500 K, 
allowing lines as weak as 0. 1 K to be detected in an hour of integration. Because the LSR velocity of each 
of the stars was known from the CO(3-2) profile, a narrower AOS spectrometer having a total bandwidth 
of 50 MHz, was used for most of these observations, providing much greater velocity resolution. The 
500 MHz spectrometer was used for R Cas, because its emission covered too large a velocity range to 
fit comfortably into the 50 MHz backend. All but one of the nineteen stars were detected in this higher 
transition. 

3. Calibration 

The receiver was sensitive to radiation in both sidebands, which were separated by 2.8 GHz. Indeed, 
when CO(3-2) was located in the upper sideband, ^^SiO(8-7) was detected in the lower sideband for 
several objects ( R Ser, R Leo and especially W Hya). The relative sensitivity of the receiver to the two 
sidebands is difficult to determine, and it changes when the receiver's mechanical tuning structures are 
moved. For this reason, the carbon star IRC+10°216 was observed whenever the receiver was mechanically 
retuned, to assess its sensitivity to the signal sideband. Corrected for the main beam efficiency as measured 



* The CSO is operated by the California Institute of Technology under funding from the National 
Science Foundation, Contract #AST-93-13929. 



-4- 



Table 1. - Survey Star Positions, Integration Times and Spectrum Noise Levels 

Nondetections are listed in parenthesis 



R.A. Decl. Time rms 

Star (1950) (1950) (min.) (niK) 



(S Scl) 


00 


12 


51 


11 


-32 


19 


22 


5 


27 


2 


89 


(T And) 


00 


19 


46 


43 


26 


43 


10 


3 


7 


1 


180 


T Cas 


00 


20 


31 


50 


55 


30 


56 





23 


3 


100 


CO(4-3) 


















54 


4 


210 


(TU And) ^ 


00 


29 


44 


70 


25 


45 


12 





13 


6 


98 


W And 


02 


14 


22 


70 


44 


04 


26 





8 


4 


75 


Get 


02 


16 


49 


04 


-03 


12 


13 


4 


5 


2 


210 


CO(4-3) 


















9 


7 


380 


T Ari ^ 


02 


45 


32 


04 


17 


18 


07 


2 


108 


7 


18 


R Hor 


02 


52 


13 


00 


-50 


05 


32 





3 


9 


300 


U Ari '^''^ 


03 


08 


16 


20 


14 


36 


40 





31 


7 


46 


S Pic ^ 


05 


09 


37 


20 


-48 


34 


00 





50 


5 


59 


R Aur 


05 


13 


15 


26 


53 


31 


56 


6 


8 


7 


110 


CO(4-3) 


















31 


1 


140 


(W Aur) ^'-^'^ 


05 


23 


31 


20 


36 


51 


40 





28 


5 


47 


S Ori ^'<^ 


05 


26 


32 


68 


-04 


43 


51 


7 


14 


2 


51 


CO(4-3) 


















69 


2 


170 


U Aur 


05 


38 


53 


50 


32 


00 


57 





12 


9 


98 


CO(4-3) 


















47 


9 


67 


U Ori 


05 


52 


51 


00 


20 


10 


06 


2 


16 


8 


94 


CO(4-3) 


















60 


8 


110 


(V Mon) 


06 


20 


11 


69 


-02 


10 


07 





19 


4 


44 


(X Gem) 


06 


43 


55 


20 


30 


19 


53 





40 


1 


25 


(Y Mon) ^ 


06 


54 


05 


80 


11 


18 


33 





50 


5 


97 


(RS Mon) ^ 


07 


04 


49 


30 


05 


03 


58 





50 


5 


150 


S CMi ^ 


07 


30 


00 


00 


08 


25 


35 





18 


1 


35 


R Cnc ^ 


08 


13 


48 


54 


11 


52 


52 


6 


14 


9 


53 


W Cnc ^'<^ 


09 


06 


57 


90 


25 


27 


06 





78 


9 


24 


(R Car) 


09 


30 


59 


21 


-62 


34 


01 


1 


29 


1 


280 


X Hya 


09 


33 


06 


92 


-14 


28 


03 


5 


45 


3 


28 


R LMi 


09 


42 


34 


40 


34 


44 


35 





10 


4 


62 


CO(4-3) 


















63 


4 


91 


R Leo 


09 


44 


52 


23 


11 


39 


41 


9 


16 


9 


50 


CO(4-3) 


















43 


2 


110 


(R UMa) 


10 


41 


07 


67 


69 


02 


22 


3 


62 


2 


30 


(X Cen) ^ 


11 


46 


41 


10 


-41 


28 


47 





34 


9 


47 


(R Crv) 


12 


17 


02 


32 


-18 


58 


41 


7 


19 


4 


79 


(R Vir) 


12 


35 


57 


67 


07 


15 


47 





16 


8 


89 


R Hya 


13 


26 


58 


40 


-23 


01 


27 





10 


4 


98 


CO(4-3) 


















20 


7 


270 


S Vir -4'^ 


13 


30 


22 


80 


-06 


56 


18 





23 


9 


47 


CO(4-3) 


















32 


4 


84 


W Hya 


13 


46 


12 


50 


-28 


07 


11 





5 


8 


160 



R.A. Decl. Time rms 



Star (1950) (1950) (min.) (mK) 



W Hya CO(4-3) 


















16 


9 


410 


(R CVn) 


13 


46 


48 


10 


39 


47 


28 





51 


8 


31 


(R Cen) 


14 


12 


56 


90 


-59 


40 


55 


2 


34 


9 


170 


(U UMi) ^ 


14 


16 


13 


40 


67 


01 


29 





16 


8 


46 


(S Boo) 


14 


21 


12 


40 


54 


02 


12 


4 


14 


2 


71 


(R Boo) 


14 


34 


59 


27 


26 


57 


09 


4 


22 





73 


(Y Lib) ^ 


15 


09 


02 


50 


-05 


49 


27 





11 


7 


76 


S CrB ^■ 


15 


19 


21 


53 


31 


32 


46 


5 


9 


1 


82 


CO(4-3) 


















63 


4 


100 


RS Lib ^'<^ 


15 


21 


24 


30 


-22 


44 


06 





24 


6 


50 


CO(4-3) 


















25 


9 


118 


(S UMi) ^ 


15 


31 


24 


10 


78 


47 


54 





19 


4 


100 


(T Nor) 


15 


40 


11 


71 


-54 


49 


44 


8 


31 


7 


160 


R Ser ^ 


15 


48 


23 


24 


15 


17 


02 


7 


24 


6 


36 


(CO(4-3)) 


















64 


7 


85 


(Z Sco) ^ 


16 


03 


04 


10 


-21 


35 


55 





6 


5 


86 


RU Her ^'<^ 


16 


08 


08 


60 


25 


11 


59 





14 


2 


79 


CO(4-3) 


















64 


7 


128 


U Her ^ 


16 


23 


34 


86 


19 


00 


18 





12 


9 


68 


CO(4-3) 


















24 


6 


110 


(S Her) 


16 


49 


37 


14 


15 


01 


27 


8 


19 


4 


68 


RS Sco ^'<^ 


16 


51 


59 


82 


-45 


01 


22 


9 


29 


1 


47 


RR Sco '^''^ 


16 


53 


26 


31 


-30 


30 


08 


3 


31 


1 


33 


(R Oph) 


17 


04 


53 


43 


-16 


01 


39 


5 


24 


6 


53 


X Oph ^ 


18 


35 


57 


51 


08 


47 


19 


7 


8 


4 


86 


R Aql 


19 


03 


57 


70 


08 


09 


08 





3 


2 


120 


CO(4-3) 


















19 


4 


300 


RT Aql ^'^'<^ 


19 


35 


39 


70 


11 


36 


24 





19 


4 


46 


X Cyg 


19 


48 


38 


30 


32 


47 


08 





6 


5 


190 


CO(4-3) 


















19 


4 


240 


(RR Sgr) ^ 


19 


52 


50 


20 


-29 


19 


21 





9 


1 


73 


Z Cyg 


20 


00 


02 


60 


49 


54 


07 





23 


3 


35 


T Cep -4 


21 


08 


53 


30 


68 


17 


12 





12 


9 


110 


(S Lac) 


22 


26 


49 


37 


40 


03 


34 


3 


13 


6 


140 


R Peg 


23 


04 


08 


23 


10 


16 


21 


5 


27 


2 


51 


W Peg ^ 


23 


17 


22 


70 


26 


00 


18 





9 


7 


80 


(S Peg) 


23 


18 


00 


95 


08 


38 


42 


1 


13 


6 


93 


(R Aqr) 


23 


41 


14 


18 


-15 


33 


41 


8 


72 


5 


36 


R Cas 


23 


55 


51 


68 


51 


06 


36 


4 


5 


2 


100 


CO(4-3) 


















15 


5 


400 


(Z Peg) 


23 


57 


33 


20 


25 


36 


29 





18 


8 


80 


(W Cet) 


23 


59 


33 


64 


-14 


57 


15 


2 


45 


3 


33 
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^Coordinates for these objects are from the IRAS PSC. 

^Four off-source positions spaced symmetrically 30 arc seconds away from the source posi- 
tion were used when these objects objects were observed, in order to minimize contamination 
from extended galactic emission, if any. 

These objects showed only weak CO emission, so they were observed and the detection 
was confirmed on at least 1 additional night. 

CVn was marginally detected when initially observed, but was not detected when re- 
observed for confirmation. 



on planets, the average value of T^b for IRC+10°216 for all tunings was 29 K. This value is in reasonably 
good agreement with the main beam temperature of 32 K, obtained by Wang et al. (1994), who employed 
a single-sideband filter in their measurements, which therefore did not suffer from any calibration errors 
introduced by unequal sideband gains. The value of Tj^b for IRC+10°216 is apt to be constant even if the 
central star's luminosity varies, since the rotational levels of CO in the envelope are populated primarily 
through collisions (Kwan & Hill 1977). For this reason, the temperature scales for each observation were 
scaled by the factor required to yield a T^b of 29.0 K for IRC+10°216. For the CO(4-3) observations, 
observations of IRC+10°216 were not available for each tuning. For these observations the antenna 
temperature was corrected only for the main beam coupUng efficiency, which was measured to be 60% 
on both Jupiter and Mars. 



4. Analysis 

Linear baselines were fitted to regions of the spectra judged to be free of emission, and were then 
subtracted from the spectra. For those objects which were detected, a model line profile was fitted to the 
data. The model (Morris 1984) assumes that CO is expelled from the central star in a spherical shell, at a 
nearly constant velocity (which is much larger than either the thermal or turbulent velocities). The line's 
temperature profile T(v) has the form 



In the above equations T^b is the fine's temperature at the systemic velocity Vlsr, and Vq is the 
outflow velocity of the CO. a is related to the optical depth of the CO transition; large values correspond 
to high optical depths and give rise to parabolic line-shapes. Small values indicate that the transition is 
optically thin, and give rise to flat-topped lines. /3 is the ratio of the size of the source to that of the 
telescope's beam. Large values of /3 give rise to profiles with "horns" at Vlsr ± Vq on optically thin 
lines, or somewhat flattened parabolae in the optically thick case. Because none of the spectra show very 
pronounced horns, and most of them are too noisy to allow discrimination between a flattened parabola 
and a perfect one, for the purposes of profile fitting all sources are assumed to be completely unresolved 
(/3 = 0). The remaining 4 parameters were adjusted for their least-squares best fit. Table 2 fists which 
stars were detected, the results of the profile analysis, and the 95% confidence limits on the derived 
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Table 2. - Line Fit Results 











VLSR 


V 

'0 


out! 




K km 


ya) 




^^Mii a ) 


T Cas 


LO ± 0.06 


17 


0.8 


-9.3 ± 0.2 


10.5 ± 0.5 


(4-3) 


0.45 ± 0.07 


5.5 


> 4 


-5.7 ± 0.6 


10.5 ±1.1 


W And 


r\ c f\ 1 f\ f\/' 

0.59 ± 0.06 


7.0 


> 4 


-36.7 ± 0.3 


8.9 ± 0.5 


o Cet 


12 ± 0.2 


75 


^ A 

> 4 


A ^ c 1 r\ r\ 

46.5 ± 0.0 


A Ci 1 r\ rv 

4.8 ± 0.0 


(4-3) ^ 


10 ± 0.15 


67.9 


1.6 


A Z' A 1 

46.4 ± 0.0 


4.7 ± 0.1 


T Ari 


0.10 ± 0.01 


0.48 


1 A 

1.4 


-0.8 ± 0.3 


2.7 ± 0.7 


R Hor 


3.1 ± 0.2 


23 


1.9 


35.7 ± 0.2 


5.2 ± 0.3 


U Ari 


0.29 ± 0.04 


1.7 


1.3 


-55.3 ± 0.2 


4.0 ± 0.3 


S Pic 


0.33 ± 0.04 


5.2 


2.5 


f\ f\ 1 rv c 

-0.9 ± 0.5 


1 1 /" 1 1 rv 

11.6 ± 1.0 


R Aur 


1.2 ± 0.07 


16 


> 4 


-3.4 ± 0.2 


9.6 ± 0.5 


(4-3) 


0.90 ± 0.04 


11.4 


> 4 


-1.0 ± 0.2 


1 rv rv 1 rv o 

10.0 ± 0.3 


S Ori 


0.45 ± 0.04 


3.4 


> 4 


14.5 ± 0.2 


5.7 ± 0.5 


(4-3) 


0.39 ± 0.07 


2.9 


2.9 


14.1 ± 0.5 


6.4 ± 0.9 


U Aur 


0.41 ± 0.04 


7.1 


0.1 


7.0 ± 0.5 


9.2 ± 0.5 


(4-3) 


0.17 ± 0.02 


2.2 


3.8 


6.8 ± 0.5 


9.8 ± 0.8 


U Ori 


1.1 ± 0.07 


11 


> 4 


n ^ 1 rv 

-37.5 ± 0.2 


7.5 ± 0.3 


(4-3) 


0.72 ± 0.04 


5.8 


> 4 


-37.9 ± 0.2 


6.5 ± 0.3 


S CMi 


0.31 ± 0.02 


2.1 


1.7 


51.7 ± 0.2 


4.7 ± 0.3 


R Cnc 


0.44 ± 0.04 


2.7 


0.3 


15.0 ± 0.2 


3.5 ± 0.3 


W Cnc 


0.097 ± 0.02 


0.73 


0.1 


o o 1 rv 

36.8 ± 0.2 


A 1 rv c 

4.2 ± 0.5 


X Hya 


0.097 ± 0.01 


0.85 


1.1 


26.2 ± 0.5 


5.9 ±1.1 


R LMi 


0.82 ± 0.02 


9.5 


^ A 

> 4 


1.3 ± 0.2 


0^7 1 r\ o 

8.7 ± 0.3 


(4-3) 


0.63 ± 0.02 


7.1 


> 4 


0.6 ± 0.2 


8.5 ± 0.2 


R Leo 


3.7 ± 0.02 


37 


> 4 


-0.4 ± 0.0 


A 1 rv rv 

7.4 ± 0.0 


(4-3) 


2.8 ± 0.04 


28.1 


3.5 


-0.4 ± 0.0 


7.4 ± 0.1 


r» TT 

R Hya 


2.5 ± 0.06 


37 


2.3 


1 r\ 1 1 f\ 

-10.1 ± 0.2 


1 r\ c 1 r\ 

10.5 ± 0.2 


(4-3) 


4.6 ±0.11 


46.1 


> 4 


-10.0 ±0.1 


8.2 ± 0.1 


S Vir 


0.42 ± 0.04 


2.5 


2.3 


9.8 ± 0.2 


4.4 ± 0.5 


(4-3) 


0.35 ± 0.04 


2.3 


> 4 


10.1 ± 0.2 


5.7 ± 0.4 


W Hya 


2.1 ± 0.07 


29 


0.3 


'^ rv rv i rv 

39.9 ± 0.2 


8.0 ± 0.2 


(4-3) 


3.1 ± 0.09 


40.3 


1.0 


41.0 ±0.1 


8.4 ± 0.2 


S CrB 


0.88 ± 0.04 


9.1 


> 4 


1.5 ± 0.2 


7.7 ± 0.3 


(4-3) 


0.60 ± 0.04 


A r\ 

4.9 


^ A 

> 4 


1.4 ± 0.2 


^ o 1 r\ o 

6.3 ± 0.3 


RS Lib 


0.35 ± 0.04 


2.4 


1.7 


7.8 ± 0.3 


4.9 ± 0.5 


(4-3) 


0.45 ± 0.06 


3.1 


^ A 

> 4 


"1 /' 1 rv o 

7.6 ± 0.3 


c c 1 rv c 

5.5 ± 0.5 


R Ser 


0.37 ± 0.02 


2.5 


> 4 


o 1 ^ 1 rv o 

31.7 ± 0.3 


5.1 ± 0.5 


(4-3) 






no? detected 






RU Her 


0.79 ± 0.04 


r\ A 

9.4 


> 4 


11/1 1 rv o 

-11.4 ± 0.3 


o rv 1 rv o 

8.9 ± 0.3 


(4-3) 


0.53 ± 0.04 


5.7 


1.1 


-11.1 ± 0.2 


7.2 ± 0.4 


TT Hp>r 

u rier 




/ . J 




1 T 7 -1- n ^ 
— vj.i It yj.j 


1 1 .J> It U.D 


(4-3) 


0.36 ± 0.02 


4.9 


0.1 


-13.7 ±0.1 


7.4 ±0.2 


RS Sco 


0.35 ± 0.04 


2.3 


> 4 


6.5 ± 0.3 


5.0 ±0.3 


RR Sco 


0.19 ±0.02 


0.85 


> 4 


-28.9 ± 0.3 


3.1 ±0.5 


X Oph 


0.33 ± 0.06 


2.9 


1.3 


-55.0 ±0.5 


6.1 ±0.9 


R Aql 


3.5 ± 0.06 


45 


> 4 


46.0 ± 0.2 


9.5 ± 0.2 


(4-3) 


3.1 ±0.09 


36.9 


> 4 


47.1 ±0.1 


9.2 ±0.2 






( continued 


on nejcf page) 
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Table 2. (continued) 





Tmb 


Area 




Vlsr 


Vo 


Star 


(K) 


Kkms"' 


(a) 


(kms-i) 


(kms-i) 


RT Aql 


0.22 ± 0.02 


1.9 


> 4 


-28.6 ± 0.5 


6.6 ±0.6 


X Cyg 


4.7 ±0.1 


63 


> 4 


7.8 ±0.2 


9.9 ± 0.2 


(4-3) 


4.0 ± 0.07 


52.7 


1.4 


9.9 ±0.1 


9.0 ±0.1 


ZCyg 


0.27 ± 0.02 


1.5 


> 4 


-147.7 ±0.2 


4.0 ±0.3 


T Cep 


1.5 ±0.07 


9.4 


> 4 


-3.0 ±0.2 


4.8 ±0.3 


R Peg 


0.23 ± 0.02 


2.5 


> 4 


22.7 ± 0.5 


8.3 ±0.8 


W Peg 


0.31 ±0.04 


3.6 


1.2 


-15.5 ±0.5 


7.7 ± 0.9 


R Cas 


4.7 ± 0.04 


87 


0.9 


23.7 ± 0.0 


12.1 ±0.2 


(4-3) 


6.0 ± 0.29 


102.4 


> 4 


25.6 ± 0.2 


13.3 ±0.3 



^ The CO(4-3) profile of o Ceti not well fit by equation 1 for any combination of parameters 
(see figure 2). 



parameters. These confidence limits were derived from the covarience matrix produced by the nonlinear 
least-squares program which was used to derive the parameters; they are merely formal errors, and do 
not include such effects as telescope pointing errors etc. The values of Vq have been corrected for the 
broadening caused by the 0.9 kms~' resolution of the 500 MHz AOS, by assuming Vq and the velocity 
resolution add in quadrature. No such correction was applied to the CO(4-3) line fits (except for that of 
R Cas), because the more narrow 50 MHz AOS used for those observations did not broaden the lines 
measurably. Figures la-d and 2a-b show the CO(3-2) and CO(4-3) spectra of the detected stars, along 
with the fitted profiles. 

Three of the stars; T Ari, R Cnc and RR Sco, have outflow velocities of less than 4 kms~\ lower 
than any Mira variable previously detected in CO (to the author's knowledge). The average outflow 
velocity is only 7.0 kms~^ This is substantially below the value found in previous surveys in which 
candidates were selected by infrared criteria. For example the survey by Nyman et al. (1992) used IRAS 
color criteria to select candidates. Their survey detected 20 oxygen-rich Miras, with an average outflow 
velocity of 13.3 kms~^ Zuckerman and Dyck (1986a, 1986b, 1989, Zuckernam, Dyck & Claussen 1986) 
surveyed the evolved stars with the largest fluxes in the AFGL and IRAS catalogs; among their detections 
were 14 oxygen-rich Miras with an average outflow velocity of 11.3 kms~^ 

4. 1 Detections and Nondetections 

IRAS PSC positions were used for 24 of the 66 survey stars. Even for stars which are bright in 
the infrared the PSC positions have a mean error of ~ 9" (Beichman et al. 1988), a sizable fraction of 
the CSO's 20" CO(3-2) beam. As can be seen in table 1, the detection rate for objects for which IRAS 
coordinates were used was not high. However the estimated distances for these objects is 50% greater 
on average than those of the stars with SAO coordinates. It is likely that the low detection rate for the 
IRAS coordinate stars is due at least in part to their greater average distance. If we consider only those 
stars more distant than 300 pc, there was not a dramatic difference between the detection rate for the 
SAO coordinate stars (46% detected) and the stars for which IRAS positions were used (44% detected). 
Nonetheless any individual non-detection must be viewed with suspicion if its position was taken from 
the PSC. 



Figure 1 □ 
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Fig. 3. The horizontal axis denotes the decimal subtype of the star's spectral type at visual 
maximum. A value of 5.5 could denote a spectral type of either M5.5 or S5.5. The bold line shows the 
percentage of stars detected. The dashed line shows the number of stars examined for each spectral 
type. 



The characteristic that best discriminates between the stars which were detected and those which 
weren't is the spectral type's decimal subtype (Yerkes system), taken from the RCMV. As figure 3 
shows, the detection percentage is a steep and nearly monotonic function of the subtype. None of the 
13 stars examined with subtypes below 5.5 were detected, while all 4 with subtypes greater than 7 were. 
An attempt was made to find an infrared color index which divided the stars cleanly into detected and 
nondetected groups. All colors of the form 

f\{i) f\(2) f\(3) fm f\(5) fm fxu) fm 

^ - fX(9) fX(lO) fXiU) fXiU) fX(n) fX(U) fX(l5) fX(l6) 

were examined, where // and fx are the 1 and K band fluxes from the TMSS catalog, and /n^ 
are from the AFGL catalog, the four long wavelength fluxes are from the IRAS PSC and < A(n) < 2. 
No color of this form was found to be able to predict which stars would be detectable in all cases. 
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Table 3. - Calculated Luminosities, Distances and Mass Loss Rates 









Dist ance 


M 


RjVIax 






Dist ance 


M 


RjVIax 


Star 


Mv 


(pc) (10- 


Meyr-') 


(10" cm) 


Star 


Mv 


(pc) 


(10-^ Meyr-') 


(10" cm) 


T Gas 


-u 




350 


110 ±36 


16. 


S Vir 


-1.87 


450 


5.6 ± 2 


5.6 


(4-3) 








45 ± 17 




(4-3) 






5.7 ± 2 




W And 


n 
-u 


1 

ID 


280 


34 ± 12 


9.7 


W Hya 


-0.12 


80 ^ 


18 ± 7 


7.5 


o Get 


-i 


Q /I 

84 


100 ^ 


36 ± 12 


14. 


(4-3) 






21 ± 8 




(4-3) 








21 ± 8 




S GrR 


-0.83 


310 ^ 


35 ± 12 


11. 


T Ari 


-i 


DO 


350 


0.6 ± 0.2 


2.4 


(4-3) 






11 ±4 




R Hor 


A 

-U 


OQ 

ao 


190 


28. ± 11 


12. 


RS Lib 


-0.11 


210 


3.7 ± 1 


3.5 


U Ari 


1 

- J 


/I /I 
44 


710 ^ 


6.5 ± 2 


6.3 


(4-3) 






3.5 ± 1 




S Pic 


-U 


1 1 


440 


59 ± 24 


11. 


R Ser 


-1.35 


330 ^ 


4.5 ± 1 


4.7 


R Aur 


-1 


i i 


450 


160 ± 60 


20. 


R.U Her 


-0.89 


410 


75 ± 24 


14. 


(4-3) 








1 10 ± 40 




(4-3) 






28 ± 10 




S Ori 


A 

-U 


i 


480 


21 ± 7 


9.4 


U Her 


-0.72 


360 


59 ± 24 


11. 










22 ± 8 




(4-3) 






14 ± 5 




U Aur 


A 

-U 


71 
i i 


500 


63 ± 24 


13. 


RS Sco 


-1.03 


300 


2.2 ± 0.8 


4.3 


(4-3) 








24 ± 9 




RR Sco 


-0.96 


200 


0.6 ± 0.2 


2.2 


U Ori 


A 

-U 


on 


240 


28 ± 10 


9.5 


X Oph 


-0.80 


270 


4.5 ± 1 


4.2 










9.1 ± 3 




R Aql 


-0.67 


190 


110 ± 36 


17. 


S GMi 


-0 


52 


340 


4.1 ± 1 


4.6 


(4-3) 






76 ± 30 




R Gnc 


-0 


88 


300 


2.3 ± 0.8 


4.0 


RT Aql 


-0.52 


450 


10 ±4 


6.2 


W Gnc 


-1 


61 


720 ^ 


3.1 ± 1 


4.3 




-0.29 


115 


82 ± 24 


14. 


X Hya 


-0 


16 


400 


3.0 ± 1 


3.6 


(4-3) 






46 ± 20 




R LMi 


-0 


47 


260 ^ 


36 ± 12 


10. 


Z Cyg 


-1.11 


490 


4.0 ± 1 


5.0 


(4-3) 








23 ± 9 




T Gep 


-0.94 


220 


9.9 ±4 


7.2 


R Leo 


-0 


46 


120 ^ 


41 ± 12 


12. 


R Peg 


-0.74 


380 ^ 


12 ± 5 


6.1 


(4-3) 








25 ± 9 




W Peg 


0.03 


270 ^ 


9.8 ±4 


5.6 


R Hya 


-0 


76 


110 ^ 


42 ± 12 


10. 


R Gas 


-0.78 


230 


290 ± 110 


24. 


(4-3) 








39 ± 20 




(4-3) 






360 ± 100 





^Distance from Jura and Kleinmann 1992, ^Distance from Gatewood 1992 
'^Distance from Stein 1991 



although several of these colors were nearly able to do so. Also, no color was a more reliable indicator 
of detectabiUty than the IRAS 15 ii flux alone. Three of the undetected stars; R Aqr, R Car and R Cen, 
were consistently included among the detected stars by all the colors and fluxes which most reUably 
discriminated between detectable and not detectable sources. R Aqr is a symbiotic variable, and this 
might explain its lack of detectable CO emission. R Cen and R Car are normal Miras however, and the 
lack of CO surrounding these stars is intriguing. Both have larger fluxes in all four IRAS bands than many 
of the stars which were detected, indicating that they do have substantial dust shells. As neither of these 
stars rise above 10 degrees from the southern horizon, it's tempting to dismiss their nondetection as a 
result of poor instrument performance at low elevations. However Nyman et al. (1992) also observed R 
Car from the southern hemisphere (at SEST), and they too did not detect it. The V band amplitudes of R 
Car and R Cen, 6.0 and 6.5 mag. respectively (Eggen 1992), are fairly low for Mira variables, but they 
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are as large or larger than several of the detected stars. However at visual band maximum the spectral 
types of R Car and R Cen, M5e and M4.5e respectively (RCMV), are earlier than any of the detected 
objects. Perhaps the hotter photospheres of these two stars prevents the dust forming near enough to 
the star to allow efficient dust-gas momentum coupling and mass loss. Alternatively, these stars may be 
losing mass at rates comparable to cooler Miras, but carbon monoxide might not be fully associated in 
their photospheres. It would be surprising if R Car and R Cen had strong winds with no CO, because 
a CO envelope has been detected surrounding a Ori (Knapp et al. 1982) which has an even warmer 
photosphere (M21ab). 

4.2 Mass Loss Rates 

After the line profile parameters had been calculated, an artificial CO(3-2) profile was calculated 
using a Large Velocity Gradient (LVG) radiative transfer program (Morris 1980). The mass loss rate was 
adjusted until the peak Tmb of the artificial profile matched that of the parameterized profile. 

The CO envelope was assumed to be truncated by photodissociation, and this radius was calculated 
using the photodissociation rate of Mamon et al. (1988). Table 3 lists the calculated absolute magnitudes, 
distances, mass loss rates (M) and the outer radii of the CO envelopes for the stars which were detected. 
The highest (R Cas: 2.9 x IQ-^M^yr-i) and lowest (T Ari and RR Sco: 6 x lO-^Mgyr"') mass loss 
rates differ by nearly 3 orders of magnitude. The star's distance, luminosity and fractional abundance of 
CO were input parameters to the LVG model, and the sensitivity of the calculated value of M to errors 
in these parameters was investigated by varying the parameters by their estimated lo- uncertainties, and 
recalculating M. 

Errors in the star's estimated luminosity and distance are not independent, because the distance was 
calculated using the luminosity provided by the period - spectral type - luminosity relationship (PSLR). If 
a particular star is actually more luminous than indicated by the PSLR, it provides more infrared radiation 
to excite CO, which leads to increased CO(3-2) radiation from a given amount of CO. Increasing the 
estimated luminosity of the star decreases the value of M calculated by the LVG model, because less 
molecular material is needed to produce the same line strength. However if the star's luminosity has 
been underestimated, then its derived distance will also be too low. Increasing the distance estimate 
given to the LVG increases the calculated value of M, because less of the CO(3-2) radiation from a 
given amount of CO will arrive at the telescope. Therefore the errors in M introduced by errors in the 
luminosity and distance will have opposite signs. The net effect is that varying the luminosity by ±0.4 
magnitudes, the mean spread in the PSLR (Cehs 1981), changes the derived distances by (+20%,-17%) 
and M by ±23%, on average. The M values shown in Table 3 were calculated with the assumption 
/ = [CO]/[H2] = 3 X 10""*, as would be the case if hydrogen and carbon were present in the solar 
abundance ratio, and half of the carbon was incorporated into CO molecules. Decreasing / to 2 x 10""* 
increases M by 39% on average, while increasing / to 4 x 10~"^ on average decreases M by 21%. If 
errors in the estimated luminosity and / are independent, then they should add in quadrature, and the 
uncertainty in the calculated value of M arising from these two sources is ± 35%. Error estimates for the 
value of M for individuals stars, using the range of parameters discussed above, are given in Table 3. 

Figure 4 shows a log-log plot of the calculated mass loss rates versus the outflow velocity. There's 
an obvious strong correlation, indicating a power-law relationship holds between these quantities. The 
best-fit power-law is 

M ^ 2.94 X 10-^0 v3-35^^_,, M^yr-i. (2) 
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Fig. 4. A Log-Log plot of mass loss rate as a function of outflow velocity is shown for the 
36 Mira stars which were detected. The line shown is the linear least-squares best fit to the data 
(equation 2). Filled circles indicate stars with integrated CO(3-2) line intensities of 3 K kms~'or 
greater, which lie above line D in figure 5. Open circles represent stars which fall below line C in 
figure 5. 



Ironically, the star which deviates most strongly from this relationship is Mira itself! The most likely 
reason for this is that Mira's true outflow velocity is substantially larger than the 4.8 kms~^ found by 
fitting equation 1 to the data. The wings of Mira's CO(3-2) profile extend to about 10 kms~^ from the 
line center (Knapp 1995). Plotting Mira in figure 4 with a higher expansion velocity would move it closer 
to the locus of points, but the peculiar high velocity wings on both the 3-2 and 4-3 profiles indicate that 
the simple spherical, constant expansion velocity LVG model used to calculate M is too simplistic for 
this star. 

It is important to recognize that selection effects could influence the appearance of figure 4 and hence 
equation 2. Figure 5 shows the integrated intensity as a function of expansion velocity, on a log-log plot. 
If we assume that the weakest lines are unresolved and optically thin, then the line area is approximately 
given by twice the product of the peak antenna temperature (T), and Vq . If the spectrometer is well 

1 /2 

behaved, the noise in the measurement of the fine area will be proportional to Vo , so the signal/noise 
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1 /2 

for a weak line will be proportional to TVq . For a fixed signal/noise detection threshold, the minimum 
detectable antenna temperature, Tmin, should therefore be proportional to Vq and the minimum 

1 /2 

detectable line area should be T^inVo . Thus one would expect the points of figure 5 to be bounded 
by a line with a slope of 1/2. However, on a more subjective level, very wide, weak lines are apt to be 
overlooked or incorrectly attributed to instrument baseline irregularities, leading to false nondetections. In 
this way one could end up with a list of detections in which all the least intense lines, which for a given 
distance indicate a low mass loss rate, have small outflow velocities. This is worrisome, because such 
a bias might lead to a relationship resembling equation 2. If all lines above some threshold temperature 
were detected regardless Vo , then one would expect the points to be bounded by a line with a slope of 
1 . If broad, weak lines were overlooked, then one would expect the bounding Une to have a slope greater 
than 1. In fact, the points all lie on or above line A in figure 5, which has a slope of 0.7, suggesting that 
weak, broad lines were not systematically rejected. However, there are no points near line A with large 
expansion velocities, and if one star (X Hya) is rejected, the points seem to be bounded more closely by 
line B, which has a slope of 1.8. So the sample is probably too small to allow the importance of selection 
effects to be quantitatively assessed. A reasonable guess would be that the 7 points falling below line C 
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Fig. 6. An IRAS 2-color diagram, for tire detected and nondetected stars is siiown. 

Tire abscissa is 2.5 logiQ(J25/i/-Fi2/i) and the ordinate is 2.51ogiQ(J60/i/J25/i). Open circles represent 
nondetections, filled circles represent detections. The size of the filled circles is proportional to the 
log of the mass loss rate listed in table 3. Five survey stars with poor quality 12, 25, or 60/tt fluxes in 
the PSC (RS Mon, R Cen, T Nor, RS Sco and R Aqr) have not been included in this plot. 



on figure 5, which has a slope of 1, represent stars with antenna temperatures which might not have been 
recognized as detections, had their outflow velocities been higher. In figure 4, these 7 stars are represented 
by open circles, and it is clear that the line produced by equation 2 is still a good fit to the data even if 

these points are rejected. Above 3 K kms~^ (line D in figure 5), there is no correlation at all between 
line area and outflow velocity. Figure 4 shows that even if only those points lying within this uncorrected 
region of figure 5 are used, there is still a correlation between M and Vq . 

Van der Veen and Habing (1988) showed that oxygen rich Miras and OH/IR stars fall along a broad 
track on an IRAS 2-color diagram constructed from the PSC 12, 25 and 60// fluxes. They further showed 
that as one moves along this path from the warmer to the cooler end, stars with increasingly large mass 
loss rates are found. These authors interpret this as an evolutionary progression. Figure 6 shows such a 
2-color diagram for the stars in this survey. Surprisingly, the detected and nondetected stars are not well 
segregated and the mass loss rate seems unrelated to either of the colors. Though not shown in figure 6, 
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Fig. 7. Figures A tiirougii D siiow tiie mass loss rate as a function of far infrared luminosity 
(calculated as flux x distance^) plotted for each of the IRAS bands. Stars with poorly determined 
IRAS fluxes are not included. In each case, the luminosity has been normalized to that of o Ceti. 
Figures E and F show there is no correlation between distance and the calculated mass loss rates 
and 60/[tm luminosity, implying that these quantities aren't dominated by the effects of errors in the 
distance estimates. 



colors involving the IRAS 100// flux were also examined and they too appeared uncorrected with the 
mass loss rate. It has been shown that for highly obscured stars the mass loss rate is strongly dependent 
on IRAS color, and empirical formulae have been constructed to give mass loss estimates based on these 
colors (van der Veen 1989). The results presented here show that such formulae are not applicable to 
stars with visible photospheres and low values of M. On the other hand, figure 7 shows that the mass 
loss rates in table 3 do correlate with all four IRAS fluxes if they are normalized to a common distance, 
although the correlation with lOOyum flux is very weak. The correlation is strongest with the 60// flux. 
Such a correlation should be viewed with caution, because errors in the stellar distance estimate influence 
both the calculated mass loss rate and the distance normalized flux in the same manor, and can produce 
a spurious correlation. However neither M nor the 60yum luminosity shows any correlation with distance 
(see figures 7E and 7F), therefore it seems likely that the correlation of M with infrared luminosity is real. 
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Fig. 8. (A) shows the distribution of CO(3-2) outflow velocities as a function of decimal 
spectral subtype at V band maximum. (B) shows the averages of the points in plot A. The error 
bars indicate the standard deviation. (C) shows the mass loss rates from table 3 as a function of 
decimal subtype, and (D) once again shows the averages. 



4.3 The Wind Velocity 



There is general agreement that the copious mass loss from heavily dust enshrouded stars, such as 
the OH/IR stars, is driven by radiation pressure on dust grains. However there is no such consensus on 
the cause of the much more modest mass loss seen in visible Miras. Netzer and Ehtzur (1993) solved 
the equations of motion for dust and gas in a radiation driven wind, and concluded that there is a lower 
limit of ~ 10~' M0yr~^ for radiation driven mass loss, well above many of the values of M in table 3. 
Fusi-Pecci and Renzini (1975) showed that acoustic energy could drive mass loss in red giant stars. Yet 
the increase in outflow velocity at later spectral types shown in figures 8A&B is what would be expected 
in a radiation driven wind. Netzer and Elitzur (1993) show that for a such a wind, if the terminal velocity 
of the gas greatly exceeds the initial velocity that the gas has at the base of the acceleration region, then 
Vo cx R,7 where Rq is the distance at which dust forms. Thus cooler stars, which should allow dust 
to form at smaller radii, will have higher values of Vq. The absence of detectable CO emission from 
the stars with spectral types earher than M5.5 indicates that either such stars have extremely low mass 



- 22 - 



loss rates or that they have mass loss rates similar to the cooler stars, but the gas is atomic. Knapp et 
al. (1995) looked for the 609yum line of C I in the envelope of Mira (spectral type M5.5) but did not 
detect it. Thus it would be surprising if large amounts of atomic carbon would be seen surrounding Miras 
with only slightly warmer photospheres. It seems more likely that these warmer Miras have very low 
mass loss rates. These warm Miras may inhibit the formation of dust throughout the region where the 
gas density is high enough to support a significant mass loss rate. It therefore appears that the results of 
this survey support radiation pressure on grains as the driving mechanism for even low mass loss rates. 
It is surprising however, that M is not correlated with spectral type (figures 8C&D). 

One of the strongest pieces of evidence that the winds from AGB stars are driven by radiation 
pressure is the fact that the momentum present in the wind is almost always less than or comparable to 
the momentum available in the radiation field (Knapp et al. 1982). Assuming radiation pressure drives 
the wind, the rate of increase of the wind's mechanical momentum (MVq) can exceed the momentum 
available from the star's radiation (L*/c) only if the light is scattered in an optically thick dust envelope. 
All of the stars in this survey must have optically thin dust envelopes, because the photospheres are still 
visible. If we assume that all Miras have a luminosity of about 3400 (the average of the oxygen 
Miras in the LMC survey of Feast et al. 1989), and if we use equation 2 to estimate the mass loss rate, 
then there is a maximum value of Vq , Vmax, for which MVq < L*/c. This value is 17 kms~^ Because 
equation 2 is such a strong function of Vo , Vmax is rather insensitive to the assumed Mira luminosity; if 
the true luminosities of Mira variables ranged from 10^ to 10^ Lq, this would only introduce a range of 
±40% in Vmax- It seems therefore that any AGB star losing mass by the same mechanism as these Miras 
will be obscured by an optically thick envelope if Vq is ^ 17 kms~^ 

In the standard radiation driven model of the wind, the dust feels the radiation pressure, and propels 
the gas outward via collisions. Thus the dust must be moving though the gas, at a velocity called the 
"drift velocity," for which Goldreich and Scoville (1976) derived the expression 



where Q is the momentum coupling constant between the radiation field and the dust. If we assume that 
the dust consists of silicate spheres with radii of 0. 1 ii, and the radiation field is that of a 2000 K black 
body, equation 3 becomes 



3ge^tL,Vo\ ^'^ /3 X 0.023 x 34OOL0 x Vo^ 



2Mc 7 V 2Mc 



where Qext, the dust extinction efficiency, has been taken from a tabulation by Draine (1987). Substituting 
equation 2 for M gives 

«d ^ 61 X V-i i* ^ ^ kms-i . (4) 

» o 

At Vo = 3 kms~^ this gives a drift velocity of 20 kms~\ which is the velocity above which significant 
ablation of dust by sputtering is believed to occur (Wickramasignhe 1972; Kwok 1975). The lowest 
expansion velocity for any of the stars in this survey is that of T Ari, 2.7 kms~'. While the assumption 
of uniform spherical grains used to derive equation 4 is a crude approximation at best, this result provides 
a possible explanation for why no Miras have outflow velocities significantly below 3 kms~^ 
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Table 4. - Comparison of Outflow Velocity Measurements for Previously Detected Stars 





CO Rotational Transition 






CO Rotational Transition 




Star 


1-0 


2-1 


3-2 


4-3 


Star 


1-0 


2-1 


3-2 


4-3 


T Cas 


10.5 ^ 


5.2^ 


10.5 


10.5 


R Leo 


4.0^ 


6.5^ 


7.4 


7.4 


W And 


11.0^ 




8.9 




R Hya 




7.2^ 


10.5 


8.2 


o Cet 


10.0 ^ 


5.0 


4.8 


4.7 


WHya 




8.8^ 


8.0 


8.4 


R Hor 


5.4 ^ 


6.6^ 


5.2 




R Aql 




10.1 ^ 


9.5 


9.2 


S Pic 


11.8 ^ 




11.6 




X Cyg 


10.2 ^ 




9.9 


9.0 


R Aur 


11.1 ^ 




9.6 


10.0 


T Cep 




5.0 


4.8 




U Ori 


5.8 « 




7.5 


6.5 


R Cas 


14.3 ^" 




12.1 


13.3 


R LMi 


6.5 ^ 


6.0^' 


8.7 


8.5 













Sources: (A) Knapp, 1986; (B) Zuckerman and Dyck, 1986a; (C) Zuckerman, Dyck and 

Claussen, 1986; (D) Nyman et al. , 1992; (E) Knapp and Morris, 1985; (F) Knapp and 
Sutin, 1989; (G) Kastner, Zuckerman, Dyck and Sopka, 1989; (H) Knapp et al. 1982; (I) 
Wannier and Sahai, 1986; (J) Zuckerman and Dyck, 1989 



Although the CO emission hnes from these stars presents a radiative transfer problem which is too 

complex to be solved analytically, it is likely that in all cases the emission from higher J transitions 
arises from a region closer to the star than does a low J transition. In the case of a large mass loss 
rate (M ^ lO~^M0yr~'), this is because the rotational levels are populated by collisions (Morris 1980, 
Schonberg 1988), which tends to force the partition function to reflect a single temperature at a given 
radial distance. This temperature will decrease with increasing distance from the star, so the Boltzman 
factor will allow high J levels to be populated only at small radh. For low mass loss rates, coUisional 
excitation is less important than radiative pumping via the CO vibrational transitions near 4.6// (Morris 
1980, Schonberg 1988). Although infrared absorption and subsequent decay can either increase or decrease 
J, the branching ratios of the transitions insure that an increase in J is more likely, therefore the net effect 
of infrared pumping is to increase the population of the higher J transitions. However the spontaneous 
decay rate of a rotational transition is proportional to J^, and as J increases the pump rate must increase 
also to maintain a significant population. This means that as J increases infrared pumping will only be 
effective at smaller radii, again implying that high J emission arises from small radial distances. Assuming 
the CO rotational ladder can be characterised by a temperature distribution of the form T{r) cx r~" , Van 
der Veen and Oiofsson (1989) derived an expression giving the approximate relative locations from which 
the different rotational emission lines arise 

Rj^j-i ^ ( 2 
iZi^o W(^ + l) 

If a = 0.7 (Kwan & Linke 1982), this imphes the CO(4-3) probes the stellar wind at a radius only ~ 1 /25 
as large as CO(l-O), and reflects the conditions in the wind only a few hundred AU from the star. 

Table 4 lists the outflow velocities, culled from the hterature, for the stars in this survey which 
had previously been detected at a lower J transition. One must be cautious in comparing the outflow 
velocities obtained by different authors who no doubt processed their data with different algorithms. Even 
for the very strong source Mira, the velocities obtained for a single transition, CO(2-l), show considerable 
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scatter ( Knapp & Morris 1985: 5.0 kms~'; Zuckerman & Dyck 1986a: 4.6 kms~h Zuckerman, Dyck & 
Claussen 1986: 8.3 kms~^; Zuckermann & Dyck 1989: 4.3 kms~'). However if the gaseous component 
of the wind is still being accelerated at radii greater than a few hundred AU, the value obtained for Vo 
should on average be lower as increasingly high transitions are examined. Table 4 shows no evidence of 
this. 

Two groups have recently studied the acceleration of SiO in oxygen-rich red giants, including x 
Cygni. Lucas et al. (1992) mapped SiO(2-l) using the IRAM interferometer, and concluded that the 
gaseous component of the envelope had attained only 1/2 of its terminal velocity at a distance of 330 AU 
from the star, and only 3/4 at a distance of 730 AU. In contrast Sahai and Bieging (1993) mapped the same 
spectral line in several of the same objects using the BIMA array, and concluded that the acceleration 
zone was only 1/10 as large as Lucas et al. had deduced. The results presented here, which show no 
systematic decrease in the the outflow velocity as higher J transitions of CO are observed, argue in favor 
of the Sahai and Bieging model. 

5. Conclusions 

A survey of stars from the RCMV within 500 pc was conducted; of the 66 stars examined 36 were 
detected. The survey data led to the following conclusions: 

1) A Mira's spectral type is the best predictor of whether an extensive CO envelope exists or not. None 
of the stars examined with spectral types earlier than M5.5 were detected, all of those later than M7 were. 

2) For the relatively low mass loss rates of these Miras there is no correlation between far infrared color 
and the mass loss rate indicated by the CO emission, in contrast to the findings of CO surveys containing 
heavily obscured stars. The far infrared luminosity is, however, correlated to the mass loss rate. 

3) Two stars, R Cen and R Car, appear to have extensive dust envelopes, yet were not detected in CO(3-2). 
Because these stars are near the southern limit of the sky observable from the CSO, it would be useful to 
examine these stars carefully from an observatory in the southern hemisphere. It would also be interesting 
to search for atomic carbon emission from these stars, to see if the lack of CO is caused by the stars' 
relatively warm photospheres, rather than an undetectablely low mass loss rate. 

4) The outflow velocity increases as stars with cooler photospheres are examined, as would be expected 
if radiation pressure on dust grains drives the stellar wind. 

5) The mass loss rate is strongly correlated to the outflow velocity. A power-law relationship holds 
between these two quantities over almost 3 orders of magnitude in M. 

6) At the low end of the observed outflow velocities (~ 3 kms~^), the dust drift velocity may be large 
enough to cause the destruction of the dust grains though dust-gas collisions. 

7) At the high end of the observed outflow velocities (~ 12 kms~^), the momentum in the wind is 
beginning to approach the total amount of momentum available in the radiation field. At velocities ^17 

kms~', the power-law Vq - M relationship indicates that all the momentum in the radiation field will be 
transferred to the wind; this can only occur if the dust envelope is optically thick in the visible and near 
infrared. 

8) An examination of the outflow velocities obtained at differing CO rotational transitions shows no 
general trend in Vq as high J transitions are observed. 
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